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ABSTRACT 

We study the effect that the dark matter background (DMB) has on the gravitational energy content 
and, in general, on the star formation efficie ncy of a molecu l ar clo ud (MC). We first analyze the effect 
that a dark matter halo, described by the INavarro et al.l (jl996[ ) density profile, has on the energy 
budget of a spherical, homogeneous, cloud located at different distances from the halo center. We 
found that MCs located in the innermost regions of a massive galaxy can feel a contraction force 
greater than their self-gravity due to the incorporation of the potential of the galaxy's dark matter 
halo. We also calculated analytically the gravitational perturbation that a MC produces over a uniform 
DMB (uniform at the scales of a MC) and how this perturbation will affect the evolution of the MC 
itself. The study shows that the star formation in a MC will be considerably enhanced if the cloud is 
located in a dense and low velocity dark matter environment. We confirm our results by measuring 
the star formation efficiency in numerical simulations of the formation and evolution of MCs within 
different DMBs. Our study indicates that there are situations where the dark matter's gravitational 
contribution to the evolution of the molecular clouds should not be neglected. 

Subject headings: ISM: clouds — evolution — dark matter 



1. INTRODUCTION 

A great effort from the astronomical community is cur- 
rently directed towards the detailed understanding of 
star formation, in particular to the precise mechanisms 
involved in the conversion of some part of a molecu- 
lar cloud's mass into stars. The evolution of molecular 
clouds (MCs) and the formation of stars within them are 
thought to be regulated by actors as diverse as gravity, 
large-scale flows, turbulence, magnetic f ields, and stellar 
feedback, among others (see reviews bv lShu et al .1119871 : 
IVazquez- Semadeni 20Tol). 

While evaluating the dynamical state of a molecular 
cloud due to all relevant forces acting on it, the Virial 
Theorem is frequently invoked to define an equilibrium 
condition by equaling the gravitational energy W of the 
cloud to twice its kinetic energy K: 



2K = -W. 



(1) 



Traditionally, this relationship has been interpreted 
in terms of clouds being in a state of quasi-equilibrium 
and long-lived entities, and deviations from this condi- 
tion are assumed to mean that the cloud is collapsing or 
expanding, depending on w hi ch term is dominant (e.g . , 
iMvers fc Goodman! Il988l bl: iMcKee fc Zwdbel [19921) . 
However, it has been demonstrated that a cloud full- 
filing the so-called equilibrium con dition will not nec- 
essar ily remain dynamically stable (Ball esteros-Paredesl 
120061 ). In fact, numerical simulations show that clouds 
collapsing in a chaotic an d hierarchical way will develop 
such a relationship (e.g., [Vazquez - Semadeni et al .1120071 : 
iBallesteros-Paredes et al .1 120 1 1 al l H) ~ 

Furthermore, when calculating the Virial Theorem for 
MCs, it is a common practice to consider that the grav- 

1 a.suarez@crya.unam.mx 



itational term is the gravitational energy, assuming that 
the cloud is, in practice, an isolated entity. However, 
different recent works have showed that the media sur- 
rounding the dense structures can importantly alter the 
evolution of such entities fe.g.. | Burkert fc Hartma"nnl 



120041: iGomez et alJl2QQ^ JVazauez -Semadeni et al.L 
iHeitsch fc Hartman 3 |2008|; [Ballesteros -Paredes et al.l 
l2009b ,a). In these works, the influence of bar ionic matter 
outside the region of interest is studied. 

On an apparently completely different topic, evidence 
from a wide range of observations has long led as- 
tronomers to argue in favor of the existence of much more 
matter than is actually seen: the elusive dark matter. Al- 
though the very nature of this mass component has not 
been discovered yet, it is presumed to smoothly perme- 
ate the majority of galaxies and in many cases to greatly 
surpass their visible extension, while being their domi- 
nant mass component. Compared to the sizes of giant 
Galactic molecular clouds, dark matter halos are huge 
and, for most practical purposes, homogeneous. Thus, 
one would expect that such a distribution will not cause 
a substantial gravitational effect on the cloud. 

In the present contribution, the idea that dark matter 
halos may have an effect on the energy content of the gas 
in galaxies is investigated. As a first approach, we focus 
our study in the center of massive galaxies, where the ef- 
fect is expected to be the greatest. A cusped dark matter 
profile such as the one proposed by INavarro et al.l (|1996l 
hereafter NFW) would provide the central region of a 
halo with a potential well which could importantly con- 
tribute to the gravitational energy of a molecular cloud 
located around there. 

In addition, we explore the idea that a barionic mass 
concentration is able to introduce a perturbation in an 
otherwise homogeneous dark matter distribution. As 
such, the mere presence of a molecular cloud would mod- 
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ify the dark matter halo in which it is embedded. This 
perturbation would in turn place the cloud in a local ex- 
ternal potential well, which would play as an extra agent 
in the struggle that determines its dynamical state. In 
this case, the external medium to the MC will contribute 
to its collapse. When compared to a situation identical to 
this one but without the dark mater background (DMB), 
one would expect that a MC in this scenario would dis- 
play an enhanced capacity to form stars, since there is an 
extra factor working in the same direction as the cloud's 
internal gravitational energy. As such, measuring the 
star formation efficiency (SFE) of clouds embedded in 
different DMBs would hint towards the importance of 
the environment in their evolution in each case. We an- 
alytically study the situation and give light as to what 
physical conditions would be required for this to be an 
important effect. 

Furthermore, we present results from numerical simu- 
lations of molecular cloud formation and evolution from 
two convergent monoatomic flows, placed in different 
dark matter contexts. Although the DMBs explored are 
not very realistic, they are designed to test the analytical 
predictions and illustrate the importance of accounting 
for the external gravitational potential in molecular cloud 
analysis. 

The structure of the article is as follows: in section 
$2] we introduce the tidal energy term that is part of 
the total gravitational energy of a matter distribution. 
In section we describe a semi-analytical procedure to 
evaluate the tidal contribution that the complete halo of 
a galaxy can have on a molecular cloud and apply it to 
a few scenarios. In section §4.11 we analytically derive 
an expression to asess the tidal effect of the dark matter 
background perturb ation caused by the molecular cloud 
itself, while section §4.21 evaluates this effect in different 
dark matter background environments. In section $5] we 
present numerical simulations of cloud formation, in or- 
der to test the analytical predictions in a more realistic 
scenario, and in $6]we show the results obtained. Finally, 
Ogives general conclusions for the two studied effects. 

2. FULL GRAVITATIONAL CONTENT OF A 
MOLECULAR CLOUD IN A GALAXY 

Molecular clouds, like any other physical system, are 
bound to follow all the internal and external forces act- 
ing on them and determining their behavior over time. 
Along the present article, the total gravitational content 
for MCs situated in different environments will be evalu- 
ated, in order to look for the effect that their surround- 
ings can have on their evolution. Hence, this section de- 
rives an expression for the tidal energy component, which 
will allow for an evaluation of its gravitational effect on 
the cloud. 

With the intention of determining a MCs dynamical 
condition, an energy balance is commonly performed to 
it by invoking the Virial Theorem. In this theorem, the 
gravitational energy term W considered for a mass distri- 
bution of volume V and density p, embedded in a gravi- 
tational potential is given by 

W = - [ pxi^- d 3 x (2) 
Jv dxi 

(see, e.g., iShu et al.ll"l987D . When applied to MCs, this 



term is traditionally assumed to be equal to the gravita- 
tional energy of the cloud: 



grav 



d 3 x, 



(3) 



where p c \ represents the density of the cloud and <£> c i 
the potential energy that it generates. This assump- 
tion considers that the relevant gravitational poten- 
tial is due only to the mass of the cloud, which is a 
valid approximation for an isolated mass distribution 
(|Ballesteros-Pared es 2006). However, in a given galaxy, 
there is a considerable amount of material (both bari- 
onic a nd dark matter) coexisting with a MC. Thus, fol- 
lowing Balle steros-Paredes et al.l (|2009bl ). one may write 
the gravitational potential $ of a MC in a galaxy as the 
potential due to the mass of the cloud itself, <3> c i, plus the 
contribution of any other external sources, ^ext, 

$ = $ cl + $ ext , (4) 

such that the gravitational term entering the Virial the- 
orem becomes 



W = E grav + Wext- 



where 



Pel Xi 



dx- 



dV, 



(5) 



(6) 



is the tidal energy, i.e., the gravitational energy con- 
tained in the cloud, but due to any mass distribution 
external to the cloud. As exp lained with detail in 
iBallesteros-Paredes et al] (2009b), the external gravita- 
tional energy can either contribute to the collapse of the 
cloud or to its disruption, depending on the concavity of 
the external gravitational potential where it is located. 

In the following sections, the tidal energy felt by a given 
MC will be compared to its internal gravitational energy; 
its relative importance will determine the role that any 
material surrounding a MC will play in its evolution and 
future star formation. 

3. TIDAL ENERGY IN THE POTENTIAL WELL OF 
A DARK MATTER HALO 

To evaluate the effect that dark matter might have on 
star formation, we first consider a toy model consistent of 
a spherical molecular cloud embedded in a spherical dark 
matter halo, with the cloud located at a distance s from 
the halo's center. The halo density profile is given by the 
simul ation- inspired NFW density profile (jNavarro et al.l 
fl996h : 



Pnfw(0 



(r/r 3 )(l + r/r s ) 2 



(7) 



where r s is a scale radius, the radius where the logarith- 
mic derivative of the density is equal to —2, and p s is a 
characteristic density, the density at r = r s . The profile 
can be expressed in a dimensionless form by 



^NFW 



(77) = rj 1 (1 +77)" 



(8) 



where ^nfw(^) = Pnfw(0/Ps an d rj = r/r s . As seen in 
Fig. [Dor eqs. and flgj), the NFW profile has a cusp 
at the center. It goes to infinity as r? -1 for 77 — > 0. 
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Fig. 1. — Dimensionless NFW density profile (<$nfw) as a func- 
tion of radius (77). The density has a cusp in the center of the 
matter distribution and diverges at 77 = 0. This profile appropri- 
ately describes dark matter halos found in numerical simulations. 

The NFW gravitational potential can be expressed as 

4>NFw( r ) = -4nG Ps r 2 s ln{1 + r/rs \ (9) 

r 

where G is the gravitational constant. A dimensionless 
version of it would be 

tfNFwfo) = (10) 

rj 

where ^nfw = — </>NFw/47rGp s r;? and rj = r/r s . It is 
noteworthy that, although the density profile is diver- 
gent, the potential is finite everywhere due to the mild 
divergence of the density; the limit of the NFW potential 
as the radius goes to zero is 

lim 0nfw -4nGp s r 2 s . (11) 

r— ^0 

We assume, for simplicity, a molecular cloud modeled 
as a sphere with constant density p c \ and radius R. When 
calculating W ex t ( e Q- EJ, it is important to note that the 
integration should be done over the volume of the molec- 
ular cloud, while the expression for the NFW potential 
is centered on the dark matter halo. A coordinate trans- 
formation is then necessary to match both systems of 
reference (separated by the distance s), which makes the 
expression significantly more elaborate for any configu- 
ration where the cloud's center doesn't coincide with the 
center of the halo (i.e. s / 0). For this reason, a nu- 
merical integrator needs to be used to obtain the tidal 
energy. 

To get a quantitative idea of the importance of the 
contribution of the dark matter halo, Wext is com- 
pared to the internal gravitational energy content of the 
cloud. For a spherical and homogeneous distribution, 
this amounts to 

E g = -^ 2 Gp 2 cl R 5 . (12) 

The energy ratio W GX t/E g measures the relative impor- 
tance of the gravitational energy contribution of the dark 
matter compared to the gravitational energy of the cloud 
itself. 



We have, in total, five parameters available to explore: 
two NFW halo parameters (characteristic density and 
radius of the density profile), two molecular cloud prop- 
erties (density and extent of the cloud), and the distance 
of the cloud from the center of the halo (s). We choose 
to fix first the halo characteristics and the distance s and 
calculate W ex _t/E g values for a range of cloud's parame- 
ters. We present energy ratio maps for typical molecular 
cloud radii and densities, one per halo type and cloud 
separation. 

The gravitational effect of a Milky Way- type dark mat- 
ter halo on a MC is first presented. As NFW parameters, 
we choose values predicted by the ACDM cosmology with 
Q m = 0.3 and ^ A = 0.7, for a halo of M vir 10 12 
(the estimated mass for our Gal axy: [Battaglia et al.l 
I2005f ). For this mass. lBullock et al.l (|2001h predict a con- 
centration c = 7.9, where c = R v { r /r s . The virial radius 
i? v i r is defined as the radius where the mean density of 
the halo is 5 times the average density of the universe. 
For the cosmology considered here S = 337 at z = 0. 
M oreover, M v \ r is simply the mass inside R v { r . The model 
bv lBullock etaTI (|200T1 ) gives the median concentration 
for a given mass, redshift and cosmology, based in sta- 
tistical results from cosmological simulations. Once we 
have M v i r and c, the NFW profile is completely spec- 
ified; for example, using the definition of R v [ r and the 
value of c one can obtain r s , while p s can simply be com- 
puted by integrating the NFW profile from to R v { r and 
using the values of r s and M v i r . Thus, for the dark mat- 
ter halo expected to permeate a Milky Way-type galaxy, 
p s = 3.6 x 1O~ 3 M pc~ 3 and r s = 25.8 kpc. 

Figure [2] shows the ratio \og(W ex _t/E g ) in gray scale, as 
a function of the cloud's density (y-axis) and radius (x- 
axis) for different values of the separation s between the 
cloud and halo's center. On the top left panel, we show 
the case where the test cloud is centered at the halo's cen- 
ter (s = pc). For this particular setup, the tidal energy 
is quite important relative to the gravitational energy of 
the cloud (the energy ratio is bigger than 0.2) for a wide 
range of cloud's densities and radii: p c \ ^ 10 3 cm -3 and 
R c i ^ 50 pc. The top right panel shows a cloud placed 
1 kpc away from the center of the halo. The importance 
of the tidal energy has already decreased for high den- 
sity clouds of the same size. As the cloud moves away 
from the center of the halo (s = 2,3 kpc, lower panels), 
the tidal contribution rapidly becomes smaller, as can be 
seen in the bottom panels. Here, the tidal contribution 
is negligible except for very low density clouds. 

Only clouds located very near the center of the halo 
feel a noticeable influence from it: the values for the 
energy ratio decay rapidly as the separation between the 
cloud and the halo increases. This is due to the central 
cuspiness of the NFW profile, shown in figure [IJ since 
Wext is dependent on the potential gradient which grows 
steeper near the origin. It is also evident that values of 
the tidal energy are more important for smaller and less 
dense clouds. Since the external gravity is "pulling" the 
particles in the cloud, it can be understood that a cloud 
of smaller diameter and lower density will be less tightly 
bound and easier to influence. 

As a second example, we present the results from an 
identical estimate but for a dark matter halo of mass 
1O 1O M . In this case, the NFW parameters employed 
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Fig. 2. — Went /Eg map for constant density clouds of radius R c \ 
and density p c \, located inside a Galactic- type dark matter halo 
(M = 1O 12 M0) and centered at a distance s from the halo center. 
When the cloud center is located at the origin of the halo, its tidal 
energy can have a comparable effect to that of the gravitational 
energy. In general, the smaller and less dense clouds feel a greater 
gravitational contribution from the external medium. 
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Fig. 3. — As in Fig. [2] but for a galaxy with a halo mass of 
1O 1O M0. A similar effect is noticed, but scaled down due to the 
potential well of such a halo being shallower. 



are p s = 8.5 x 1O _3 M pc -3 and r s = 3.8 kpc. As in 
Fig. [2j in Fig. Owe present the log(W ex t/^) maps 
for this halo's characteristics. We notice that for a halo 
with these parameter values the effect is even weaker, and 
fades away faster as the molecular cloud moves away from 
the center of the halo. This smaller contribution can be 
explained as a result of the shallower potential well of 
this less massive halo. 

It is important to mention that the ratios presented 
here only show the importance of the dark matter tidal 
energy compared to the gravitational energy of a cloud. 
We do not make any comparison with the gravitational 
influence of other possible matter components, such as 
a galaxy's central black hole or star population. With 
the results shown, we conclude that a dark matter halo 
is only capable of influencing the evolution of small, dif- 
fuse clouds which are located near the center of massive 
galaxies. 



4. GRAVITATIONAL FEEDBACK FROM A 
PERTURBED HALO 

4.1. Perturbation Analysis 

To continue with another perspective of the dark mat- 
ter influence in a molecular cloud's evolution, this and 
subsequent sections study the gravitational feedback ef- 
fect that a perturbation on the local dark matter back- 
ground (DMB), produced by a barionic cloud, has upon 
the energy budget of the cloud itself. 

Dark matter interacts with barionic matter only grav- 
itationally; as such, dark matter particles behave as a 
non-colissional fluid that feels a gravitational attraction 
to a mass distribution like a molecular cloud. Because 
galaxies are immersed in dark halos, a MC residing in it 
will coexist with a DMB. Following a dark halo distribu- 
tion like that of NFW, this DMB would have an almost 
constant density at scales of the MC, because the varia- 
tion scale of the halo's density profile is much larger than 
the extent of even a giant molecular cloud (GMC). How- 
ever, the mere presence of the barionic cloud perturbs 
the otherwise homogeneous DMB and promotes a den- 
sity enhancement in it, which generates a potential well 
that will influence the MC in return. We analyze how im- 
portant this perturbation can be to alter the evolution of 
a molecular cloud in different dark matter environments. 

Although m olecular clouds are highly irregular 
and s tructured (jComb es 1991: iBallesteros-Paredes et al.l 
120071 and references therein), we model them as Plum- 
mer spheres in order to follow up an analytical expression 
which may be easy to use. The mass distribution of such 
a sphere follows a density profile given by 



Pel 



3M cl 

47T<2 3 



-5/2 



(13) 



where M c \ is the total mass of the cloud, r is the radial 
distance from the center of the cloud, and a is the scale- 
radius of the Plummer potential, which gives an idea of 
the size of the core of the distribution. The associated 
gravitational potential can be written as 



$d(r) 



GM C 



1 



a y/l + r 2 /a 2 ' 



(14) 



where G is the gravitational constant. A cloud with these 
profiles is then placed within a DMB with constant den- 
sity pdm an d a velocity disper si on <td m- 

Following iHernandez fc Leel (12008). the dark matter 
particles are described by a Maxwell-Boltzmann velocity 
distribution function fo(v) oc exp(— v 2 /2a^ M ). For the 
purpose of this study, the dark matter particles will have 
an average velocity of 0, that is, the DMB as a bulk is 
considered at rest with respect to the barionic matter. 
Although the gravitational potential of the dark matter 
halo $o is constant at the scales of the molecular cloud, 
the MC induces a perturbation <E>i, so that the total po- 
tential $ is the addition of both contributions, i.e. 



$(r) = $oO) +e$i(r), 



(15) 



and the original density po is enhanced by a perturbation 
pi and should become 
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p(r) =Po + epi(r). 



(16) 



Similarily, the original Maxwell-Boltzmann distribu- 
tion function fo(v) of the local dark matter halo gets 
perturbed by an amount f\ (r, v) and the new distribu- 
tion function will be given by 



f(r,v) = f (v) +efi(r,v). 



(17) 



In order to know how the dark matter potential will be 
affected by the presence of the molecular cloud, we must 
solve the Boltzmann equation for non-collisional systems 



df 
dt 



-£ + v V/- V$- V„/ = 0, 



(18) 



in which V represents the gradient operator with respect 
to spatial coordinates and V v the gradient operator with 
respect to velocity coordinates. The quantities involved 
in eq. ([15]) are the ones taking into account the pertur- 
bation caused by the barionic matter, i.e., the enhanced 
potential (eq. [TSj) and distribution function (eq. [T7J) . For 
simplicity, we look for a first order stationary solution to 
the Boltzmann equation (i.e. df /dt = 0). By using the 
Jean's swindle (<3>o = 0), the fact that fo does not explic- 
itly depend on the spatial coordinates (dfo/dxi = 0), 
and considering only first order terms, the radial solu- 
tion can be obtained as 



vV/ = V*-V„/, 

= V($ + *i)-V„(/o + /i) ) 
= V$i-V„/ + V$i-V„/i, 
«V*i-V„/ . (19) 

We can calculate the gradient of the gravitational po- 
tential perturbation $1, which is actually the cloud po- 
tential $ c i, as 



ggi _ GM cl 
dr 



(l + r 2 /a 2 ) 3 / 2 ' 



(20) 



Using the explicit form of the Maxwell-Boltzmann veloc- 
ity distribution function to get the velocity gradient, 



dfo v 

-7— = 2—f0{V), 



(21) 



and after inserting the last two expressions in eq. (flU|) . 
it becomes 



rGM cl 



dfi 

' dr a 3 (l + r 2 /a 2 ) 3 / 2 



'DM 



fo(v), (22) 



which furthermore, can be integrated over velocity space 
as 



^ = -GM cl 
or 



a 3( 1+r 2/ a 2)3/2 I ^ 



p (v). (23) 



Remembering that the original DMB density has a con- 
stant value of /?dm, this equation can be solved to yield 



Pi 



GM C ; 



'DM 



^l + r 2 /a' 



=Pdm, 



(24) 



which is the density perturbation in the dark matter dis- 
tribution produced by the molecular cloud. Such a den- 
sity perturbation in the dark matter halo will produce 
its own non-constant gravitational field, which may, in 
principle, influence the gravitational energy budget of 
the molecular cloud. Our interest lies, thus, in the grav- 
itational feedback of this enhancement over the cloud 
itself. 

In order to calculate the tidal energy produced by the 
perturbed dark matter halo over the cloud, it is necessary 
to calculate the gradient of the gravitational potential 
0dm associated to the perturbation pi given by eq. ([24]) . 
The gravitational potential as a function of the radius r 
that is produced by the perturbation p\ will be given by 



pi(r')r'dr' >, 



(25) 



integrated over the radial coordinate r' and where the 
second integral will be evaluated up to some radius i?DM, 
such that the perturbation is fully contained within this 
radius, i.e. Rbm is taken such that p\(r = Rbm) = po- 
Introducing ([24]) and integrating, we get an expression 
for the perturbed dark matter potential: 



. 2nG 2 M cl a 
pDM{r) 2 Pdm 



{ (r/a)- 1 In (r/a + y/l + (r/a) 2 ) 

DM I V J 

y/l + (r/a) 2 + 2y/l + (i? DM /a) 2 ), (26) 



whose gradient is given by 

2 M cY a / a ln ( r / a + V( r / a ) 2 + x ) 



2 PDM* 



dr 



DM 



2r 2 



2r 



(27) 



Introducing this equation and the Plummer density dis- 
tribution (|13p into eq. we obtain 



27rG 2 aM 2 / 9 DM / 2R/a 



ln (r/. 



DM 



1 + (R/a) 2 

WW) 



arctan (R/a) 



(1 + (R/a) 2 ) 3 / 2 



(28) 



where R is the truncation radius of the Plummer profile 
for the molecular cloud. This expression quantifies the 
tidal energy of the perturbed dark matter distribution 
over the cloud; in order to evaluate its importance in 
the process, it should be compared to the gravitational 
energy of the cloud itself, given by equation ([3]). Us- 
ing eqs. ([T4]) and ([T3]) to evaluate this expression, the 
gravitational energy of the cloud becomes: 
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-Sgrav — 



3GM 2 
2a 



R/a 



R/a 



15 [ 



8(i + (w) A(i + (R/ a yy 



arctan (R/a) >. 



(29) 



Thus, ratio of the gravitational to the tidal energy is 



We: 



j g v 10pcy \M pc~ 



where we have defined 

f{r/a) = -1.8025 

with 



0"dm 
kms _] 



fext(R/a) 
fg(R/a) ' 



/(A/a), 
(30) 

(31) 




1 



20 30 
R/a 



40 



50 



Fig. 4. — Ratio / = f e *t/fg for a Plummer cloud of truncation 
size R and characteristic scale a = 10 pc. Note that the ratio is 
always positive, implying thus that the tida l energy always con- 
tributes to the collapse, according to eq. (1301 ) 



/extC^/a) = 



2R/a 



+ (R/a) 2 
(R/a + v 7 ! 



arctan (R/a) 



Wa)'< 



(1 + (i?/a) 2 ) 3 / 2 



(32) 



and 



7 (ii/a) 



i?/a 



i?/a 



8 (1 + (R/a) 2 ) 4 (1 + (i?/a) 2 
1 



arctan (R/a) j. (33) 

Equation ([30]) tells how important is a given dark mat- 
ter environment on the gravitational energy budget of a 
molecular cloud. Several points must be stressed from 
this equation. First of all, the ratio of tidal to gravita- 
tional energy does not depend on the mass of the cloud, 
neither on its density. Furthermore, the tidal energy on 
the cloud produced by the perturbed halo will be more 
important for larger clouds, and larger background halo 
densities. On the other hand, for a halo's larger velocity 
dispersion odm? the perturbation must be smaller and, 
as a consequence, the tidal energy will also be smaller, 
when compared to the gravitational energy, as showed 
by the cr^L dependence in eq. ([30]) . 

In Fig. 00 we plot the ratio / vs R/a, as given by eqs. 
([3TH33]) . Since the function is always positive, we note 
that the ratio of tidal to gravitational energy, eq. ([30]) . 
will also be positive for any physically consistent com- 
bination of parameters R/a, Pdm, c"dm- This means 
that the tidal energy caused by such a DMB will al- 
ways contribute to the collapse of the MC perturbing 
it. This effect, when relevant, is expected to be reflected 
on the evolution of the cloud. In a molecular region form- 
ing stars, a rise in the compressional forces acting on it 
would show an increased star formation efficiency (SFE), 
as more stars will form faster when there is relatively less 
support against collapse. In the following section we ex- 
plore the parameter space relevant to this scenario (R/a, 
Pdm, ctdm). 



4.2. Semi- analytical Results 

In order to get an idea of the magnitude of the ra- 
tio W e xt/Eg in different dark matter environments, we 
need to provide dark matter backgrounds to the molec- 
ular clouds under study. As in section §3, we make use 
of the NFW model to estimate the structural parameters 
of the dark matter halos, where the MCs will be embed- 
ded. We need to determine the halo's density ^dm and 
velocity dispersion <tdm, given its total mass M v i r and 
the radius r at which these properties are evaluated, at 
a specified redshift z. We only consider here relatively 
nearby galaxies; that is, we assume redshift z = 0. 

To find pdm and <tdm we thus proceed as follows: for 
a given halo mass M v j r , we obtain the c oncentration pa- 
rameter c following iBullock et al.l (j200l[ ). With these val- 
ues, we calculate its density profile parameters as spec- 
ified in section Once specified, we can evaluate the 
profile at the desired radius r to know the local density 
Pdm- On the other hand, in order to obtain <tdm , we 
adopt the relation given by lLokas fc Mamon I (j200lL see 
their eq. 14) for the case where velocities of dark mat- 
ter particles are considered isotropic. Such expression 
is dependent only on c v i r , M v i r and R v [ r and so the ve- 
locity dispersion obtained for a given halo is a constant. 
These recipes are used in our analysis to get the DMB 
parameters needed in different cases. 

To begin with, we first choose parameters consistent 
with a Milky- Way type dark matter halo evaluated at 
a radius of 8 kpc, adequate to the solar neighborhood. 
However, in this case we preferred to use a more realistic 
concentration parameter c = 12 . closer to the observa - 
tionally estimated value by, e.g., Battaglia et al.l ([20051 ). 
With these conditions, the procedure just described pre- 
dicts a DMB with density pDM ~ O.OlM pc -3 and a 
velocity dispersion (Tdm = 165 kins" 1 . With these pa- 
rameters, equation ([30]) predicts that in order to get a 
non- negligible tidal effect (W ext /E g ~ 0.15) the cloud's 
scale radius a has to be huge (a = 1900 pc). A cloud 
of this size greatly surpasses the extent of local molec- 
ular clouds. On the other hand, since the energy ratio 
Wext/^grav varies as a 2 , we note that for more typical 
clouds, a ~ 10 pc and R/a ~ 10, the contribution of the 
tidal energy is negligible. 
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Fig. 5. — W e xt/E g ratio for a cloud of Plummer radius a = 10 pc 
and different truncation radii (R), embedded in a dark matter halo 
typical of dwarf galaxies (M ~ 1O 7 M0), at a radius r = 0.1r s . In 
this case, even small molecular clouds would feel a a non-zero tidal 
effect from their environment. 



Figure [5] presents a more interesting case in which the 
energy ratio is non-negligible for a range of cloud sizes 
R/a. The cloud has a Plummer radius of a = 10 pc and 
it is embedded in a dark matter halo with mass M = 
1O 7 M at a distance of r = 0.1r s . We can see, in this 
case, that even small clouds of the order of 10 pc will feel 
a non-zero tidal effect from the dark matter environment. 

Figure [6] shows, in grey scale, the energy ratios 
W ex t/Eg evaluated using a NFW dark matter density 
profile with a range of velocity dispersions <jdm and den- 
sities Pdm- In this figure, the parameters of the Plummer 
cloud are a =10 pc, and R/a = 10. Similar maps are pre- 
sented in Figs. [71 and [8j but for clouds with R/a = 5 and 
1, respectively. As a reference, in these images we also 
indicate the typical densities and velocity dispersions of 
dark matter halos of galaxies with total masses ranging 
from 10 T to 10 11 M . Density profiles are computed at 
r = 0.1r s . 

Two points can be highlighted from Fig. [6] First, the 
tidal energy of the dark matter halo becomes important 
for large clouds (R ~ 100 pc) in galaxies with masses 
ranging from 10 7 to 10 8 M . Secondly and as noticed 
previously, since the ratio W ex t/E graiV scales as a 2 , this 
contribution become important also in larger galaxies for 
large enough clouds. For instance, as can be inferred 
from Fig. [7J ~ 10% of the total gravitational energy con- 
tent of a giant cloud complex (a ~ 100 pc, R ~ 0.5 kpc) 
embedded in a 10 10 M dark matter halo will be due to 
the tidal compression energy produced by the potential 
well induced by the cloud itself. 

The high similarity between Figs. [6] and [71 is likely due 
to the mass of a Plummer sphere being concentrated in 
the inner 5 characteristic radii (93% of the mass of the 
distribution lies at r < 4a). On the other hand, Figure 
[5J which shows the case R/a = 1, does show smaller 
contributions from the environment to the energy budget 
of the cloud than those shown in previous figures. The 
reason can be traced again to the inherent concentration 
of a Plummer sphere: as we truncate the sphere at r = a, 
we are only considering the inner core of the cloud, which 
has a higher gravitational energy and in which the tidal 
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Fig. 6. — W ex t/E g map for Plummer clouds with characteristic 
radius a = 10 pc and R/a=10, embedded in a homogeneous dark 
matter background of parameters pdm and (Jdm- The tidal energy 
can contribute to the collapse of clouds immersed in dwarf galaxies. 
As the energy ratio increases as a 2 , also massive galaxies can show 
this effect for large enough MCs. 
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Fig. 7. — Same as in Fig. [6] but for R/a=5. The differences with 
Fig. [6] are not significant because Plummer spheres contain most 
of their mass within 5 characteristic radii. 



contribution has more trouble to act. 

To further study the "big cloud" scenario, a different 
exploration of parameters is shown in Figure [9j in 
which the density and velocity dispersion of the dark 
matter background are set according to the theoretical 
values of a Milky Way- type halo, at a radius of 50 kpc 
(Pdm = 0.0002 M©pc~ 3 , (jdm = 150 kms -1 ), while the 
cloud characteristic radius a and truncation radius R 
are varied. This allows us to identify the characteristics 
of a mass distribution that would suffer a relevant 
contraction by the tidal forces: the figure shows that, for 
a radius larger than ~ 8000 pc (a = 8000 pc, R/a = 1, 
or a = 6000 pc, R/a = 1.3) and located at some 50 kpc 
from the center of the Milky Way, it will produce a 
perturbation in the dark matter halo which will in turn 
have a small contribution to the gravitational energy 
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Fig. 8. — Same as in Fig. [6] but for R/a=l. In this case, the 
tidal effect is diminished when compared to Fig. [6] because we are 
considering a more concentrated sphere in which the tidal forces 
have more difficulties to act. 
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Fig. 10. — Same as in Fig. but for a Milky Way dark 
matter halo at 2 kpc from its center (pdm = 0.04 Mqpc -3 , 
(7dm = 150 kms -1 ). In this case, a large cloud complex would 
suffer a tidal effect from the dark matter, compared to the previ- 
ous case. 




6000 

a[pc] 

Fig. 9. — W e xt/E g map for Plummer clouds of characteristic 
radius a and truncation radius R, embedded in a homogeneous 
dark matter background with density and velocity dispersion cor- 
responding to a Milky Way dark matter halo at 50 kpc from its 
center (pdm — 0.0002 Mqpc -3 , <t dm = 150 kms -1 ). A big mass 
distribution (like a satellite galaxy) would feel a small tidal energy 
contribution. 



it contains (around 1%). This location and dimensions 
are similar to that of a satellite galaxy like the Large 
Magellanic Cloud. The effect is indeed low, but the 
theoretical implications it suggests are interesting, as 
we could say that at least some contribution to the star 
formation seen in orbiting galaxies or in galaxy mergers 
could come from the tidal energy caused by a dark 
matter perturbation. 

Still, one more scenario to consider is that of a cloud 
closer to the center of a galaxy Milky- Way type. In Fig- 
ure [lOj the dark matter parameters are fixed to fit the 
Milky Way's dark matter halo at 2 kpc from its center 
(Pdm = 0.04 M pc -3 , ctdm = 150 kms -1 ), while the 
cloud dimensions are explored. Note that in this case, a 



large cloud complex of a ~ 1 kpc and R ~ 2 pc will have 
a tidal energy (~ 6%) that will slightly contribute to its 
collapse. This possibility suggests that a complex like the 
molecular ring in the Mily Way may be affected somehow 
by the tidal energy produced by the perturbation in the 
halo that the cloud complex imprints on it. 

5. NUMERICAL SIMULATIONS OF CLOUD 
FORMATION WITHIN A DM HALO 

Different studies support the idea that molecular 
clouds in the Milky Way are transient entities produced 
by compressions fr om large scale streams in the i n- 
terstellar medium ([B allesteros-Paredes et al. 1999a.b; 
Hartmann et al.l l2QQll ; [V azquez- Semade ni et al.l 120071: 
Heitsch fc Hartmannl I2QQ8I: I Vazquez- Semadeni et al.l 



On the other hand, in sections 14.11 and 14.21 we 
derived semi-analytical calculations which suggest that 
a dark matter well produced by the perturbation of 
a barionic distribution of mass can significantly add 
to its gravitational budget, and thus contribute to its 
collapse, if the right conditions are given. In order to 
test this scenario within a more realistic setting, we have 
performed a series of numerical experiments in which 
a dense cloud is formed by the collision of two streams 
of diffuse gas. The colliding gas cools down, becomes 
denser and collapses. While doing so, the collapsing 
cloud perturbs the locally homogeneous dark matter 
halcB in which it is embedded, producing a potential 
well. This will give an extra tidal energy contribution 
to the original cloud, which is expected to translate into 
a larger star formation activity. Given the differences 
between the semi-analytical and the numerical circum- 
stances (initial conditions and treatment), we do not 
expect semi-analytical results to faithfully reflect the 
behavior of the gas and dark matter in the simulations; 
still, the gravitational feedback effect must appear in 

2 As the scale radii of the halos considered in these simulations 
are much greater than the size of the clouds, this is a good aprox- 
imation. 
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TABLE 1 

Variable parameters in the simulations 



Pdm ctdm 
(M pc- 3 ) (kms- 1 ) 



log(M vil 
[M ]) 



r/r s 



Linestyle 
in Fig. [12] 





0.1 
0.1 
0.1 
0.1 
0.17 
0.24 





37 
10 
2 
1 

3.16 
3.16 



(...)- 
10 

8 

(•••) 6 

(•••) b 
7 
7 



(...)« 
0.15 
0.26 

0.13 
0.1 



solid 
dashed 
long-dashed 
dash-dotted 
dotted 
dot-long dashed 
short-log dashed 



Note. — Dark matter density (pdm) and vel. dispersion 
(<jdm) are shown, as well as the mass of a galaxy (M v i r ) whose 
characteristic NFW halo represents these values at the radius 
(r) given, when applicable. 
a Reference run 

b There are no realistic NFW dark matter halos with this char- 
acteristics 

the simulations. This effect can be then estimated in 
the numerical experiments through the comparison of 
the cloud's star formation efficiency with and without 
dark matter. 

To perform the simulations, we use the Adaptiv e Re- 
finement Tree (jKravtsov et al.lll997UKravtsovl2Q03l ) code 
(ART), which employs an TV— body algorithm to solve 
the gravitational interactions of the system. The chosen 
version of the code additionaly solves gas hydrodynamics 
and can han dle radiative cooling and s tar fo rmation, as 
presented bv I Vazquez- Semadeni et all (j2010f ). 

The initial conditions for the gas are taken 
from previous star form ation simulations by 
I Vazquez- Semadeni et al.l (|2010f ). We consider a 256 pc- 
per-side box filled with tenuous gas (no = 1 cm -3 ) 
at a temperature of To = 5000 K, which represents 
a warm neutral hydrogen medium (molecular weight 
/io = 1.27). The gas has an inital turbulent velocity 
fluctuation distribution of magnitude v rms = 0.1 kms -1 
everywhere. On top of that, two gas streams moving 
towards each other along the x— axis are created by 
adding an initial velocity component of 7.5 kms -1 , to 
the material contained in two 112 pc long cylinders with 
radii of 64 pc, lying on each side of the box. Since the 
sound speed for the gas amounts to c s = 7.4 kms -1 , 
the flows are trans-sonic. The gas that comprises the 
streams constitutes less than 20% of the total mass in 
the box. The cloud evolution is followed during 40 Myr. 

The local dark matter background is represented by 
particles that interact only gravitationally, both with 
other dark matter particles and with the gas. Table [T] 
shows in the first and second column the values of the 
density pdm and the velocity dispersion <jdm, respec- 
tively, used in each simulation. These were calculated 
following the schemes described in section [6] In the third 
column, we show the typical mass of a halo (galaxy) as- 
sociated with those pdm and <jdm values, given at the 
radius shown in the fourth column. 

To decouple the effect of the dark matter back- 
ground on the behavior of the gas from other effects, 
we need to study its evolution without a dark mat- 
ter background. This r eferen ce run was taken from 
I Vazquez- Semadeni et al.l (j2010h (SAF0, which stands for 
Small Amplitude Fluctuations without stellar feedback) 



and it has the same values of the gas parameters de- 
scribed above. These authors report the formation of a 
"central cloud" of irregular shape with a radius of the 
order of 10 pc, where star formation is concentrated. Al- 
though the physical configuration of the simulation dif- 
fers from the one of the analytical case we studied pre- 
viously, we believe it is a good test case for investigating 
the general effect: the promotion of a higher star forma- 
tion rate. 

In order to quantify the differencees in the evolu- 
tion of the different runs, we calculated the Star For- 
mation Efficiency (SFE) as d efined by, for example, 
I Vazquez- Semadeni et al.l (j2010[ ): 



SFE 



M* + Mr 



(34) 



where M* represents the mass contained in stars and 
M gas accounts for the mass in the dense gas (n ^ 
10 2 cm -3 ), at any given moment during the evolution 
of the system. We also compute the relative difference 
between the SFE of any run with dark matter and the 
reference run (without it), at any given time, A; that is, 



SFE re f — SFE r 



SFE 



rcf 



(35) 



In this context, the first point to clarify is whether 
these variables allow us to distinguish between the evo- 
lution of runs with different physics (e.g., no dark matter 
included vs. a particular configurations of dark matter 
halos), and the evolution of runs with the same physics, 
but different initial fluctuations. Thus, in Fig. [TTJ we 
show the time evolution of the SFE (upper panel) for 
three runs with identical physics (with dark matter prop- 
erties pdm = O.24M0pc -3 , (Jdm = 3.16 kms -1 ), but 
with different phases for the initial density and velocity 
fluctuations. However, in this case and unlike explained 
before, what is drawn as A in the lower panel of Fig.fTTlis 
the relative difference between the SFE of the run plotted 
with a solid lin43 and the runs plotted with a dashed or 
a dotted line (which differ only in the phase of the initial 
conditions). As can be seen in this figure, the differences 
among the three runs are not large, so the different initial 
conditions do not constitute a significant source for varia- 
tions in the SFE. Thus, if any of our performed runs with 
dark matter produces substantially larger differences, it 
must be due mainly to the gravitational influence of the 
dark matter. 

6. NUMERICAL RESULTS 

As commented above, the collision of the streams in the 
simulations develops instabilities in the atomic gas which 
changes its phase and becomes molecular material with 
lower temperatures and higher densities; this new phase 
forms condensations that grow with time. As the shocked 
region increases its mass, the cloud starts contracting. 
The star formation starts at ~ 25 Myr after the begining 
of the simulation, and has a rapid growth soon after. 
Around 10 Myr later, most of the gas in the simulation 
box has been transformed into stars and hence the SFE 
approaches its maximum possible value of 1 (see Fig. [T2]) . 

3 The seed used to run this model is the same used to run the 
models with dark matter shown in Tabled 
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Fig. 11. — SFE and differential values (with respect to the solid 
line run) for three different runs with different random phases of the 
initial density and velocity field, but otherwise identical physical 
conditions: pdm = O.24M0pc - ,<7dm = 3.16 km s -1 . As can be 
noticed, the differences are not large. Thus, for runs with dark 
matter halo, any difference substantially larger than these shown in 
this figure will be considered as representative of the gravitational 
influence of the dark matter halo. 

In Figure [12] we show the time evolution of the SFE 
(upper panels) and A (lower panels) for the different runs 
outlined in Table [TJ as compared to the reference run, 
which is shown as a solid line in all panels. On the left 
side we present the runs with density /?dm = O.lM pc -3 
and velocity dispersions <jdm = 37kms _1 (dashed line) 
and (Jdm = lOkms -1 (long-dashed line). Aside from the 
initial fluctuations, these runs show no significant differ- 
ence in the evolution of the SFE, as compared to the ref- 
erence run (solid line). This behavior is consistent with 
the results predicted for W ex t/ E grav by our semi-anlytical 
approach (see Figs. [6j[7]and[8]). As far as the SFE is con- 
cerned, these simulations proceeded as if the dark mat- 
ter was not present. The central panel depicts two more 
runs with the same density (/?dm = O.lM pc -3 ) but with 
lower velocity dispersions: (Jdm = 2kms _1 (dash-dotted 
line) and <tdm = lkms -1 (dotted line). In this cases, 
there is a marked difference in the overall evolution of 
the SFE. The <jdm = 2kms _1 run distinguishes from 
the previous ones in that its SFE is consistenly higher 
than the one of the reference; there is a period of time 
in which it has an SFE a factor of two higher than the 
reference value. The dotted line in the same panel has 
an SFE that is systematicaly higher than all previous 
runs. It can clearly be seen that runs with lower dark 
matter velocity dispersions (density remaining constant) 
produce higher SFEs. As can be deduced from Figure [7J 
an important contribution is expected from the external 
gravitational potential towards the collapse of the cloud: 
a ratio W ex t/E g in the range 0.3-0.8 for these cases. This 
translates into a significant enhancement of the star for- 



mation rate. The right panel, in the same Fig. [12j shows 
two other cases. These runs have both a higher density 
than those runs considered in the left and middle pan- 
els. One run has p = O.17M pc -3 (dot-long dashed line) 
and the other one p = O.24M0pc -3 (short-long dashed 
line). They are simulated with the same velocity disper- 
sion (Jdm = 3.16 kins" 1 . These models also show an 
enhanced SFE in comparison to the reference run, but 
here the contribution due to the DMB seems somewhat 
larger than the expected from the analytical predictions 
(see Figs. [6j [71 and [8]). Such a deviation is, however, 
not surprising: the predicted external contribution comes 
from an ideal model that considers a spherical molecular 
cloud at rest, while the simulation contains more realistic 
conditions. 

7. DISCUSSION AND CONCLUSIONS 

We investigated the effect that the gravitational contri- 
bution from dark matter, in the form of tidal energy, has 
on the energy balance of a molecular cloud and its evo- 
lution. On considering the large-scale dark matter halos, 
where galaxies are supposed to be immersed, we noticed 
that there is an important gravitational contribution only 
for clouds very near to the center of these halos, spe- 
cially for low-density, small molecular clouds. The effect 
of this collapse-promoting force dissapears very rapidly 
as the cloud location is moved away from the center. To 
our knowledge, this effect has not been studied before, 
perhaps in part because clouds were not expected to be 
affected by the presence of the dark matter. In light of 
these results, these predictions can be applied to struc- 
tures encountered in the central region of galaxies. In 
fact, iGenzel et al.l (j 19901 ) report on a series of observa- 
tions regarding the rich gas structure in the inner 10 pc 
of the Milky Way, comprised of a couple giant molecular 
clouds, a circumnuclear disk and an HII region. These 
entities lie at the right position, are small and clumpy 
(averaging — 3 pc in diameter) and span a larg e range 
of densities (jVollmer. B. fc Duschl. WTII 12002ft . Such 
gas distributions can, according to our results, be feel- 
ing a great gravitational compression from the external 
medium, relative to their self-gravity. These results also 
suggest that the intense star formation activity found in 
the innermost regions of large galaxies (like ours) may 
have some contribution from tidal compression caused 
by the dark matter halo of the galaxy. 

In this particular case, the presence of the central 
black hole (with an e stimated mass Mbh = 4 x 1O 6 M0; 
iSchoedel et al. 1 12002) is only relevant when the molecu- 
lar cloud structure reaches distances smaller than ~ 1 pc, 
in accordance with the radius of influence of a central 
black hole = GMbh/^o (wh ere vp is the velocity 
dispersion of stars at the center; iPeebles II1972L for our 
Galaxy, « 0.5 pc). Following the same analysis pre- 
sented here to set an example, we estimate that the gravi- 
tational energy contribution (Wrh) that the Milky Way's 
central black hole excerts on a 3 pc diameter cloud with 
density n = 10 3 cm -3 , located 5pc away from the center 
of the Galaxy is given by Wbh/E 9 < 0.01, which is not 
significant in comparison to the dark matter contribu- 
tion. 

On the other hand, the local dark matter background 
permeating a molecular cloud is shown to be able to 
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Fig. 12. — SFE and A values for the simulation runs described in Tabled The reference run is displayed as the solid line in all three 
panels. The left panel compares high velocity dispersion DMBs (cdm — 10 kms -1 ) with the reference run. Both runs in this panel have 
a density of pr>M = O.IMqpc -3 : the dashed line shows the run with (7dm = 37 kms -1 , while the long-dashed line depicts the one with 
(7dm — 10 kms -1 . No significant difference is found from the reference run, as expected from the analytical predictions. The central panel 
shows two more runs with the same density (pdm — O.IMqpc - 3 ) but smaller velocity dispersions: (7dm = 2kms~ 1 (dash-dotted line) and 
(7dm — 1 kms - 1 (dotted line). These cases show a systematically higher SFE than the reference run, also as expected. Two more runs with 
higher densities are displayed on the right panel: pr>M = O.17M0pc -3 , (7dm = 3.16 kms -1 (dot-long dashed line) and pdm = O.24M0pc -3 , 
(7dm — 3.16 kms -1 (short-long dashed line). Again, the star formation seems a lot more active than the reference run. 



contribute to its star formation if the appropriate con- 
ditions are fulfilled. The high dark matter densities and 
low velocity dispersions, found in dwarf spheroidal galax- 
ies, needed for the tidal energy to be important for the 
evolution of molecular clouds are not expected in galax- 
ies where there is a significant amount of molecular gas; 
namely, spiral and irregular galaxies. Yet, small galax- 
ies dominated by dark matter with a significant fraction 
of gas are not ruled out. Moreover, although the Solar 
Neighborhood does not meet the dark matter conditions 
where star formation would be greatly enhanced by its 
presence, we predict a small contribution for big clouds. 



This work has shown that despite the huge difference 
between the length scales of molecular clouds and dark 
matter halos, the latter could influence the evolution of 
a star forming region. We show that this effect is negli- 
gible in our local environment, but that it can have pos- 
sible applications in other environments, other locations 
or other galaxies. In summay, our study firmly suggests 
that the external dark matter gravitational contribution 
should be considered, when analyzing the dynamical evo- 
lution of a molecular cloud, if the right conditions are 
fullfilled. 
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